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Abstract. We study the variational convergence of a family of two-dimensional Ginzburg- 
Landau functionals arising in the study of superfluidity or thin-film superconductivity, as the 
Ginzburg-Landau parameter e tends to 0. In this regime and for large enough applied rotations 
(for superfiuids) or magnetic fields (for superconductors), the minimizers acquire quantized 
point singularities (vortices) . We focus on situations in which an unbounded number of vortices 
accumulate along a prescribed Jordan curve or a simple arc in the domain. This is known to 
occur in a circular annulus under uniform rotation, or in a simply connected domain with an 
appropriately chosen rotational vector field. We prove that, suitably normalized, the energy 
functionals F-converge to a classical energy from potential theory. Applied to global minimizers, 
our results describe the limiting distribution of vortices along the curve in terms of Green 
equilibrium measures. 
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1. Introduction 

The Ginzburg-Landau theories have had an enormous influence on both physics and 
mathematics. Physicists employ Ginzburg-Landau models in modeling superconductiv- 
ity, superfluidity, and, more recently, for rotating Bose-Einstein condensates (BECs), all 
systems which present quantized defects commonly known as vortices. In mathematics, 
starting with the work by Bethuel, Brezis & Helein [B] , many powerful methods have been 
developed to study the physical London limit, i.e., as the Ginzburg-Landau parameter 
e tends to 0. This limit corresponds to the Thomas- Fermi regime in BEC, and to an 
analogous regime in superfluids where the characteristic length scale e is very small. In 
a two-dimensional setting, vortices are essentially characterized as isolated zeroes of the 
order parameter carrying a winding number, and in the London limit as point defects 
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where energy concentration occurs. The question of whether energy minimizers develop 
vortices, where they appear in the domain, and how many there should be (for given 
boundary conditions, constant applied fields or angular velocities) has been analyzed in 
many contexts and parameter regimes. 

In this paper, we will focus on the following Ginzburg-Landau energy, arising for 
instance in the physical context of a rotating superfluid. Considering a bounded simply 
connected domain I? C M^, a smooth vector field F : ^ M'^, > and £ > 0, we 
define the functional 

u e H\V-C) ^ T,{u) ^ {^iVi^P + ^(1 - - • dx . 

Identifying with the complex plane C, we denote by 

i{u) ■.= u/\\/ue L^{V;R^), 

the pre-Jacobian of u. The L^-vector field j{u) may also be written as j{u) = (m, Vu), 
where (•,•) is the standard inner product of two complex numbers, viewed as vectors 
in 

In the case of uniform rotation, that is V{x) = x'^ = {—X2,xi) and with D a 
disk, Serfaty [T7] studied minimizers of a closely related functional (see Remark II. 5p 
to determine the critical value fii = Sli(e) of the angular speed at which vortices 
first appear (see also |10lll| for BECs). She finds that minimizers acquire vorticity at 
fii = k{'D)\ lne|+0(ln | lne|) for an explicitly determined constant k{'D). In a series of pa- 
pers, culminating with the publication of the research monograph [16 , Sandier & Serfaty 
developed powerful tools to study vortices in Ginzburg-Landau models. Although they 
primarily work with the full Ginzburg-Landau model with magnetic field, the methods 
apply as well to the functional J-^ above. In particular, their results apply to the near- 
critical regime in simply connected domains. In our setting, their results show that for any 
simply connected domain V, the first order expansion of the critical value Qi for vortex 
existence in minimizing configurations is also given by fc(P)| lne| for some constant k{T>). 
Moreover the locus of concentration of vortices for = ili + o(| lne|) is given by the set 
of maxima of with C, the solution of the following boundary- value problem: 

f-AC = curlV^ inP, 

[C ondV. ^ ' ' 

The constant fc(I?) is then determined by 
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2|C|max 

where |C|max denotes the maximum value of If, for instance, V is real-analytic and 
curlV^ is nonnegative, then so is the solution C, and the maximum is gcnerically attained 
at a finite number of points in T) (see e.g. [7|). In this situation, if 17 = fJi -|- o(| lne|), 
minimizers exhibit vortex concentration at isolated points, and the number of vortices 
remains uniformly bounded whenever il — J7i is of order 0(ln | Inej), see |16I17I10I11) . 

The case of a multiply connected domain provides a slightly different qualitative 
picture. In a work on rotating Bose-Einstein condensates, Aftalion, Alama & Bronsard [1] 
considered a similar functional in a domain given by a circular annulus A (centered at the 
origin) and again with uniform rotation V{x) — x^ (see Remark ll.51) . Unlike the simply 
connected case, minimizers in the annulus may have vorticity without vortices, as the 
hole acquires positive winding at bounded rotation fJ. Then point vortices are nucleated 
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inside the interior of ^ at a critical value fii, again of leading order |lne|. Solving 
equation (jl.ip in the annulus A, one finds out that the set of maxima of the function C 
is given by a circle inside A (see Example 15. ip . Hence one can expect that, rather than 
accumulating at isolated points, vortices concentrate along this circle in the limit e — s- 0. 
The main feature proved in [T] is that if ~ fii + 0(ln | lne|), then vortices are indeed 
essentially supported by a circle S and that the total degree of these vortices is of order 
In I Inej. In other words, in the limit e — s- 0, infinitely many vortices concentrate on S, a 
phenomenon that we call vortex concentration along a curve. However the question of the 
distribution of the limiting vorticity around the circle was left open. Subsequent results 
of Alama & Bronsard |2|3j extend the result of [1] to multiply connected domains and to 
the full Ginzburg-Landau model with magnetic field and pinning potential. In contrast 
with the previous case, concentration on curves might not be a generic phenomenon for 
the Ginzburg-Landau model with magnetic field in a general multiply connected domain. 
Indeed, in this setting the vector field V represents the electromagnetic potential and it 
is an unknown of the problem. The results in |2|3j show that, for near-critical external 
applied fields, the concentration set of vortices is also given by the set of maxima of a 
certain potential related to V. This set may contain finitely many points and/or closed 
loops. Assuming that it contains closed loops, they prove that vortex concentration along 
a curve occurs, but the determination of the limiting vorticity was again left open. 

To effectively separate the question of the nature of the concentration set from the 
question of localizing vortices, we instead start with a simply connected domain T), and 
we prescribe the function ^ with ( > in V and (\dv = 0, in such a way that C is 
maximized on a prescribed curve S CC V. Then, we define 



as our vector field. We will prove that vortices will be forced to accumulate on S as e 
tends to 0. The curve E can be either a smooth Jordan curve or a smooth embedded 
simple arc, compactly contained in V. In this setting, we shall resolve the problem of 
distribution of vortices along curves, both for niinimizers and in the more general setting 
of F-convergence. In the last section we will show that in a multiply connected domain, 
and for more general vector fields V, the problem does not differ too much in nature, 
and that a similar analysis can be performed. 

To state our main result we must give more specific hypotheses on and the angular 
speed ft. We assume that ( satisfies the following assumptions: 

(HI) C e Cl^'^iV), C > in V, and (max := max^^^ Ci^) > 0; 

(H2) Ti :— {x ^ V : C,{x) — Cmax} CC P is a Jordan curve or a simple embedded arc 
of class . 

We further assume that — n{e) is near to the critical value needed for the presence of 
vortices. More precisely, 



for some function u : (0, -t-oo) (0, -foo) satisfiying uj{e) — > -l-oo with | lne| ^uj{e) — > 
as £ ^ 0. 

For u G H^{'D; C) we consider the rescaled functional 




|ln£| 



(1.2) 
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and for a nonnegative Radon measure /i on I?, we define 

^ J JVy.V 

where the function G denotes the Dirichlet Green's function of the domain P, i.e., for 
every y G G(-, y) is the solution of 

|G(-,?/) = ondV. 



(1.3) 



Our main result deals with the F-convergence of the family of functionals {F^}^-^q as 
e — > 0, and it is stated (as usual) in terms of the vorticity distribution given by the weak 
Jacobian, that is half the distributional curl of the pre-Jacobian (see e.g. [16]). 

Theorem 1.1. Assume that (HI), (H2) and ()1.2|) hold. Let En 0+ he an arbitrary 
sequence. Then, 

(i) for any {un\nen C H^{V; C) satisfying sup„ F^^{un) < +oo, there exist a subse- 
quence (not relabelled) and a nonnegative Radon measure /x in H~^(T)) supported 
by E such that 

1 .... ... .^0.1. 



■curlj(w„) > ^ strongly in {C„' (V))* ; (1.4) 

(ii) for any C H^(T>;C) such that p.4p holds for some nonnegative Radon 

measure fi in H~^{T)) supported by T,, we have 

liminf iv„(M„) > I{p) - CmaxM(2?) ; 

n — ^+oo 

(Hi) for any nonnegative Radon measure fi in H^^^D) supported by there exists a 
sequence {unjneN C H^{'D;C) such that (|1.4[) holds and 

lim (u„) /(^) - (max ^J'{^^) ■ 

n — ^+00 

As it is well known, the F-convergence theory is well suited to study asymptotics in 
minimization problems (see e.g. p ). In this context, we shall derive from Theorem ll.il 
the following convergence result for the vorticity of global minimizers, and hence solving 
the problem on the limiting distribution of vortices along E, see Remark 11.31 below. 



Corollary 1.1. Assume that (HI), (H2) and (jl.2p hold. Let e„ 0"*" be an arbitrary 
sequence. For every integer n G N, let Un € 7J^(I?;C) be a minimizer of F^^. Then, 



■curlj(u„) — > -TTT l^* strongly in {Cq' {V))* , 



Ci;(e„) n^+oc- 2/* 

where fi^ is the unique minimizer of L over all probability measures supported on E, and 
/* := /(m*)- 

Remark 1.1. As a direct application of the results in [16] (see also [H]), we shall 
see in Section 2 that for configurations {we} with i^^-energy uniformly bounded from 
above, the vorticity distribution curlj'(w£) can be approximated (with respect to the 
(Cq'"'^ (2?))*-topology) by a measure of the form 27r J^i^i diSa^ for some finite set of points 
{oijig/^ C T> and integers {di\i^i^ C Z. In other words, each point Oi can be viewed as 
an approximate vortex with winding number di. Thus the integer Z?^ = X^ie/ M«l may 
be refered to as to approximate total vorticity of the configuration u^. It is commonly 
known that approximate vortices carry a kinetic energy essentially greater than or equal 
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to ttD^\ lne| (see Section 2 for more details). With such an estimate in hand, and using 
the arguments of Section 3, we actuaUy obtain a more refined lower bound for the energy 
than the one given by Theorem II. li claim (ii). More precisely, one has 

liminf (^j^ \\^u,\^ + ^(1 - \u,\^f dx - lne|^ > , 

and 

liminf I V^C • j{ue) dx + T:D,\\ne\] > -(max ^T^) ■ 

As a consequence, if {u^} is any recovery sequence (in the sense of (in) of Theorem ll.ip . 
the liminf 's above become limits, and equality holds in each case. In analogy with [5], 
we may then say that /(/i) plays the role of renormalized energy. 

Remark 1.2 (Minimizers). From CoroUarv 11.11 and Remark 1 1. 11 we deduce that if 
is energy minimizing, then 

D,^^iu{e)+o{iu{e)) as e ^ , 

and from Theorem I l.li the minimal value of the energy expands as 

min u;^{£)F,^-^u;^ie)+oiLj\e)). 

Remark 1.3 (Equilibrium measures). The value /(/x) gives the electrostatic energy 
of a positive charge distribution fi on the set E CC V. The minimizer of / over all 
probability measures on S is called the Green equilibrium measure in associated to the 
set S, and gives the equilibrium charge distribution of a charged conductor inside of a 
neutral conducting shell, represented by dV. The value 1//* is refered to as to the capacity 
of the condenser (S, 91?). The interested reader can find in [14| many results on the 
existence and general (regularity) properties of the equilibrium measures as well as some 
examples. For instance, if I? is a disc and E is a concentric circle, then the equilibrium 
measure fi^. is the normalized arclength measure on E, see |14 [ Example II. 5. 13], and thus 
vortices are asymptotically equidistributed along E as e — > 0. However for an arbitrary 
curve E, the distribution is of course non-uniform in general. In case where E is an 
embedded arc, it is even singular at the endpoints, see [IH Example II. 5. 14]. 

Remark 1.4. In the present results the structure and regularity assumptions on the 
set E given in (H2) are mainly motivated by the physical context of [l|2|3j . However 
it will be clear that (H2) can be relaxed into weaker statments. More precisely, the 
proof of Theorem 1 1 . 1 1 relies on (H2) only for the F — limsup inequality, i.e., claim (Hi). 
The construction of the recovery sequence (see Section could be applied with minor 
modifications if the set E is for instance a finite union of piecewise arcs/ Jordan curves. 
Actually E could even have a more general structure such as a non-empty interior. In 
this later case we assume that 9E is made by finitely many arcs and Jordan curves of 
class C^. Then, given a nonnegative Radon measure /i G H^{T)) supported by E, one can 
construct a recovery sequence for /i applying the approximation techniques of Section [4] 
to /xL9E and a more standard regularization procedure for /iLint(E) as in [16] . 

In Section [5] we will show how to apply the method to a multiply connected domain 
and a more general vector field V . For simplicity, we shall consider only domains A which 
are topological annuli, i.e., A = 'D\B, where V, B are simply connected and B CC T). For 
multiply connected domains and/or a general field V, there is an extra step involved in the 
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analysis. Indeed, the leading order term in the minimal energy (of order | Inep) is due to 
the curl-free part of V which induces a diverging phase in any minimal configuration, and 
to the vorticity in the hole B which acts as a sort of giant vortex (by analogy with [1 , ) . In 
other words, to get information on the internal vorticity, one has to perform a second order 
F-convergence analysis. As in [2] we first describe minimizing vortexless configurations 
(i.e., energy minimizers over §^-valued maps) which nonetheless have vorticity around 
the hole. Then we show that for arbitrary configurations, the energy due to the curl-free 
part of V and the hole decouples nearly exactly (see Proposition 15. 2p . After separating 
out this contribution to the energy, the residual energy functional resembles above, 
and the F-convergence analysis can be done in much the same way, with some care taken 
to control the residual vorticity around the hole. The result is stated in Theorem 15.21 
The F-limit again involves a Green energy for measures supported by a prescribed set. 

Remark 1.5 (Related funtionals). As mentioned earlier, the methods here may also 
be applied to other Ginzburg-Landau functionals which have the same structure as J-^- 
A simple variant is 




This energy has been studied in more than one context. Serfaty considers in |il7j the 
minimization of this energy for rotating superfluids, with V — x-^ and under a Dirichlet 
boundary condition u\dx> = 0- The minimization of the same energy under natural (Neu- 
mann) boundary conditions also arises in a simplified model of thin-film superconductors 
introduced by Chapman, Du & Gunzburger [8] (see also Alama, Bronsard & Galvao- 
Sousa [1]). In this setting, ngV{x) = A^{x) represents the magnetic vector potential of 
the externally applied magnetic field h^x — i^eCnrW. With fi^ satisfying ()1.2|) . and a 
sequence {u^} in the energy regime of Theorem ll.il fwhich holds for minimizers under a 
homogeneous Neumann boundary condition), the two energies agree very closely, 

Te{ue) ^ Teiu,) + ^ \V\^dx + o{l) ase-^0. 

In the case of the homogeneous Dirichlet boundary condition some care must be taken 
since a singular boundary layer arises near dV as £ — > (see [12]), but otherwise the 
result of Theorem 11.11 should remain essentially the same. 

For annular domains, another functional which resembles J-^ has been used in the 
modeling of rotating BECs in certain anharmonic traps (see [I]), 



Ge{u) = 




Here V{x) — x'^ is the velocity field of uniform rotation, and the function a{x), positive 
in A, gives the trapping potential which contains the condensate. It is shown in [I] that 
^g-minimizers develop as e — + 0, infinitely many vortices concentrating along a circle 
for sufficently high rotation fi^. We believe that similar results to the ones in Section [5] 
should hold, but the method is not directly applicable here. Indeed, our analysis is based 
on "global energy estimates" of [TB] (see (12. ip in Proposition 12. ip and the presence of 
the inhomogeneity a{x) requires local estimates. Moreover the analysis of vortices for 
the energy is complicated by the fact that a{x) vanishes on the boundary, and some 
delicate estimates are required so as not to lose too much information near the boundary. 

One may also consider the full Ginzburg-Landau model of superconductivity for com- 
plex order parameter u : T> ^ C and magnetic vector potential A : "D ^ M.'^. For the 
Ginzburg-Landau model in a simply connected domain with constant applied magnetic 



T -convergence of 2D Ginzburg- Landau Junctionals 7 



field, Sandier & Serfaty [121 Theorem 9.1] have proven a F-convergence theorem of the 
form of Theorem 1 1.1 1 for applied fields of the form h^x — Hqi + (^{s), where in that case 
it is appropriate to take ln|lne| ^ uj{e) <C |lne|. This problem exhibits vortex concen- 
tration at points, and the limiting energy is obtained by rescaling around the points of 
concentration. The rescaled vorticity measures of minimizers converge to an equilibrium 
measure associated to a different problem in potential theory, a "Gauss variation" prob- 
lem whereby the charges are to be optimally placed in subject to an applied electric 
field (see [S]). Concentration on curves is possible in multiply connected regions (see 
[2|3] ). The methods described here should apply in this situationl once adapted to the 
magnetic setting, although the energy obtained as a F-limit should be for the Helmholtz 
(and not the Laplace) Green's function, —AxH{x,y)-\-H{x,y) — Sy{x). However, for the 
full Ginzburg-Landau functional with magnetic field in general domains, it is an interest- 
ing open problem in PDE to determine which (if any) non-symmetric multiply connected 
domains exhibit vortex concentration on curves. 

We finally mention a recent paper by Kashmar |13j exhibiting concentration on a circle 
for a magnetic Ginzburg-Landau functional in the disc with an inhomogeneity a{x) as 
above described by a radial step function (modelling for instance a superconducting body 
made by two different species). Here again we believe that similar results should hold but 
the method does not directly apply due to the inhomogeneity of the Ginzburg-Landau 
energy density. 

The plan of the paper is as follows. In Section 2 we prove assertion (i) of Theorem ll.il 
Section 3 tackles part (ii) oi the theorem. The upper bound of statement (Hi) is derived in 
Section 4, completing the proof of Theorem ll.il The proof of CoroUarv ll.ll is presented at 
the end of Section 4. Section [5] sketches how the preceeding arguments must be modified 
to treat annular domains. 

Notations. For any open set B gT) and any admissible map u, we denote by 



Ee{u,B) := 




the so-called Ginzburg-Landau energy of u in _B. Given e > we define 

:= {x e P : dist (x, &D) > e} , 
and for a sequence £„ —> 0+, we shall write tOn '■— w(eri) and f2„ := f^e„. 

2. Compactness of normalized v^reak Jacobians 

This section is devoted to the proof of claim (i) in Thcorcm ll.il The key ingredient to 
prove compactness of normalized weak Jacobians is the so-called "vortex balls construc- 
tion" taken from [TBI Theorem 4.1]. 

Proposition 2.1. There exists a constant Eq > such that for any < s < Eq and any 
u e H^{V;C) n C^iV) satisfying E^{u,V) < ~Je, the following holds. For any CqE^^^ < 
r < 1 there exists a finite collection of disjoint closed balls Br = {B(ai, pi)}i^i^ such 
that, writing Bi := B{ai,pi), 

(i) r = ^Pi; 

(ii) {|1 - > e^'^]r\V. C X>e n y B, =:V ; 
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(Hi) setting di :— deg(u^ dBi) if Bi C V^, and di := otherwise, 



Eeiu, V) > irDr { In f— ] - Ci ] , (2.1) 



r 



where Dj, := \di\ is assumed to be positive; 

(iv) Dr < C2\\ne\-'^E,{u,V) ; 

and Co, Ci, C2 are universal constants. Moreover, i/Co£^^^ < ri < r2 < 1 and Bi, B2 
are the corresponding families of balls, then every ball of Bi is included in a ball of B2 ■ 

In the remainder of this section, we consider an arbitrary sequence e„ — > 0+. In 
the following lemma, we prove an upper bound on the Ginzburg-Landau energy for a 
sequence having an -energy uniformly bounded from above. It will allow us to apply 
the previous proposition to such a sequence. 

Lemma 2.1. Assume that (HI) and (jl.2p hold. Let {u„}„gN be a sequence in H^(T>;C) 
such that sup„i^e„(u„) < +00. Then there exists a constant C independent of n such 
that ||'u„||^4(-x)) < C and E^^^{un,'D) < C|ln£„p. 

Proof. Observe that 

00 / |w,i||VM„|(ia: + 0(a;„) 
Jv Jv 

<\j^ \^Ur.\Hx + IIVCIIL ^ \Un\^dx + 0{loI) . (2.2) 

In particular, 

' (|u„|4 - 2(1 + 24172 |,v(||^)|^„|2 + i)dx < Oielcol) , 

so that ||u„||i4(x)) < C for a constant C independent of n. Inserting this estimate in 
I), the announced result follows easily. □ 



The first step in proving compactness of the normalized Jacobians is to show that 
the approximate total vorticity is bounded by the excess rotation ojn- We emphasize that 
here S could be any compact subset of V. 

Proposition 2.2. Assume that (HI) and p.2p hold. Let {wnjneN be a sequence in 
H^(T);C) n C^(I?) such that sup„ -F^^ (m„) < +00. Then there exist rn —> 0"*" and a 
sequence of families of balls Br„ = {B{af , pi^n)}ieir„ o,s in Proposition \K7i such that 
Dr„ ^ Cujn for some constant C > independent of n. 

Proof. Let r„ := |lne„|~''. In view of Lemma [2.11 we can apply Proposition 12 . 1 1 to m„ 
with r — r„ and n large enough. For each such n we denote by B^^ — {i3(a", pi^„)}ig/^^ 
the corresponding family of balls. For convenience we write /„ :— Ir„, Bf := B{af, pi^n), 
di_n ■= deg(u,9i?") if C di^n — otherwise, and £>„ :— D^.^. 

Step 1. We first claim that 



/ V-^C • j("n) dx = -27rrj„ d^ 



,C(ar)+o(l) asn^+00. (2.3) 
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Indeed, applying IF, Theorem 6.1] with the balls {B^}i^i^ we derive the estimate 



cm-1 j(u„) - 27r ^ di^n 



< Cr„(l + S,„K,P)) <C\\iiSn\-\ (2.4) 



thanks to Lemma O Then (g^l) follows since ri„ = 0(| lne„|) and C G (^"'^(P). 

Step 2. Without loss of generality we may assume that LOn ~ 0{Dn), otherwise there 
is nothing to prove. In view of claim (iv) in Proposition 12.11 and Lemma 12.11 we have 
Dn < 0(| ln£„|). Next, combining the lower bound (|2.ip with (|2.3p . we infer that 



0{u;l) > LolF.Sun) > vrA^(| lne„| - Cln | lne„|) - 2^ri„ ^ d,,„CK) + 

+ / |Vw„pdx + o(l) (2.5) 



as n — > +00. 

Let us now fix a sequence r/„ — s- 0+ such that 

max{ijj„, In I lne„|} = o(77„| lne„|) as n ^ +oo . (2-6) 

We group the vortex balls into the following classes and we define: 

-D* := ^ dj,„ , /,* {i G /„ : di,„ > and C(a") > (max - Vn} ] 

-^n = X! ' {j e /„ : dj,„ > and ({a^) < Cmax - ?7«} ; 

Observe that = £>* + £>+ + D~ . We claim that 

<(j max{u;„,ln|ln£„|}D„ 
" ~ 77„|lne„| 

n- / niax{^„,ln|ln£„|}j:>„ 

"""-^ u^^t;^ ' ^^■^> 

Dn<C max{t^„ , In | In e„ | } , (2.9) 

for a constant C > independent of n. In particular, if cj^ = o(|lne„|), then we can 
choose rin satisfying in addition max{u;^, (In | lne„|)^} = o(ry„ | In e„ | ) , and consequently 
yields -D+ = = for n large enough. 



We evaluate the lower bound for each class of vortex ball separately. First, we use the 
explicit form of ri„ (see (11.21) ) and the bound (^{x) < Cmax to obtain, 

7:D*J ln£„| - 2T:n„ ^ d.^^XK) > ~2TTu;„Cm..D*„ . (2.10) 

For negative degrees we have the simple estimate 

TTD-\\nen\-2TTnn ^ d^^nCK) > ^^5" | In £„ | . (2.11) 
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Then for the vortex baUs staying away from E, we have 

TT^+I ln£„| - 2^f7„ '^^.nCK) > (^1 lne«| - 27rr!„Cmax) D+ + 

-27rcj„Cniax + 27rfJ„r7„) > 

>CQ.aVnD+, (2.12) 

since w„ = o(?7„fJ„). We now insert (|2?T0| . ((2?TT|) and (|2TT2| into ((2?5)) . 



O(cj^) > ^I?„ (I ln£„| - Ci ln| ln£„|) - 27ra« V d^.nCiaT) + [ \Vu,,\^dx + o(l) 

(2.13) 

-r I I V u,„| uj. . 

Rearranging all terms we derive 



> -nCiDn In I In e„ I - 27rw„Cmax-D;; + nD^ | lne„| + Crjn^nD^ + 

+ / |Vu„pda;. 

JV\OB^ 



I ln£„| + D+7jA < C {uJnDl + In | ln£„| Ai + ^l) 
< C max{[j„ , In I In £„ I }£)„ , 

which proves ^1} and (EH). 

To prove (|2.9p . we argue as in [51 pg. 58-60] to obtain a constant C" > (independent 
of n) such that 

[ |Vu„pdx > C"L»2 . (2.14) 

Accepting (|2.14p wc return to the lower bound (|2.13p to deduce 

Dl - C (iUn + In I In £„ I) A. < 0(^1) , 

so that Dn < Cmax{a;„, In | ln£„|} and estimate (|2.9p is established. 

It remains to show (|2.14p . To this aim we identify an annular band lying outside of 
E and use the fact that the total degree is approximately constant in that band. Since 
the boundary dT> is assumed to be smooth, there exists < So < ^dist (E, dV) such that 
the function g{x) := dist (x, dV) is smooth in V^"^ with 

-.^{xeV : £„ < g{x) < So} ■ 

The level sets 

Ct:^{xeV: q{x) = t} , 

are smooth and diffeomorphic to dV for all t £ [0, So]- Define the set r„ C [0, So] by 

Tn {t e (£„, So] : Ct n D.ei^B^ - 0} . 

By the choice of r„ and claim (i) in Proposition 12.11 Tn is a finite union of disjoint 
intervals and £^((£„, So] \ Tn) < 0{\ ln£„|^^). From claim (ii) in Proposition [O] we can 
define the degree of m„ on Ct for every t £ Tn, i.e., 

Dn{t):=deg(^,Ct)^:^ f ^AX7r(^)dH\ 

\\Un\ ) 2tT Jc^ \Un\ \\Un\ J 
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where Vr denotes the tangential derivative along Ct oriented counterclockwise. Setting 
In{t) '■= {i ^ In '■ > Pi,n + 1} {pi,n being the radius of the ball B^), we have 

ieI,^{t) 

Using (|2.7p and (|2.8p we infer that for n large enough, 

Dn{t) > D„ - 2D- - D+ > ^D,, for every t e T„ . 

Denote w„ := u„/|u„|. In view of claim (ii) in Proposition l2.1[ using the Coarea Formula 
and Jensen Inequality we can estimate for n large enough, 



/ \VUn\''dx > [ |u„P|Vf„|2dx > \ r' ( I \yVn\^dn'\dt > 

Jv\UB^ Jvil\uB^ ^ Je„ \JCt / 

which completes the proof of (|2.14p . 

Step 3. If In I In £„ I < o{un) the conclusion follows from l[2Jl) . ((2?8l) and ([2Jl) . If w„ < 
0(ln I lne„|) we must refine our lower bound by growing the vortex balls. First observe 
that in this regime, (|2.7p and (|2.8p ensures that 



Dn = = (2.15) 

for n large so that each ball C I?e„ carries a nonnegative degree and D„ = D* . 

We choose a new radius s„ :— e~v^ and thus Sn > r„ for n large enough. We now 
reapply Proposition 12.11 with r = s„ to obtain a new family of larger balls {^"}jGj„, 
each new ball -B" containing one or more of the smaller balls {Sf }. By claim (ii) in 
Proposition 12.11 and (|2.15p we have Dg^ — Dn- Using the lower bound (|2.ip together 
with (|2.3p and (|2.15p . we can argue as in Step 2 to derive 

0{ojI) > ujIF.M > ^Dn (^In (t^) - Ci) - 27rn„CmaxA. + C'Dl 

> TtD,, (^In (^-^^ - Ci^ - 27rc^„Cmaxi^n + C'-D.^ . (2.16) 

Next we distinguish two cases. First assume that lnl?„ < 0(a;„). In this case, (|2.16p yields 
the inequality D^ — CuJnDn < O(w^) (with C > independent of n) so that Z?„ < 0(w„) 
as claimed. If a;„ = o(lnZ?„), we obtain the bound D^ — CD^ lnZ)„ < 0{u!n) which also 
yields Dn < 0{LUn), and the proof of Proposition 12 . 21 is complete. □ 

We are now ready to prove claim (i) in Theorem 11.11 

Theorem 2.1. Assume that (HI) and ()1.2p hold. Let En 0+ and let {un}neN be a 
sequence in H^{T>; C) such that sup„ F^^ (m„) < +oo. Then there exist a subsequence (not 
relabelled) and a nonnegative Radon measure /i supported by S such that 

Hn := — curlj(u,i) > /i strongly in (C?.'^(T>)Y . 

OJn n— >+oo 

Proof. Step 1. We first assume that {u„}„gN C H^{T>] C) D C^{'D). Using the notations 
of the previous proof, we apply Proposition 12.21 to obtain the family of vortex balls 
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{Bflig/^. Define the measure 



Since D„ < 0(cj„) we have |/2„|(I?) < C for a constant C independent of n. Therefore, 
up to a subsequence, /i„ ^ as n ^ +cx) weakly* in the sense of measures on T) for some 
finite Radon measure [i. We claim that /i is nonnegative and supported by E. Indeed, 
decompose /!„ in its Hahn decomposition, i.e., write /!„ = /2+ — /I~ where /i+ and — /x~ 
are respectively the positive and the negative parts of Then we have 

, 27rD,7 
/i„ (I?) — > as n — > +00 , 

thanks to (|2.8p . and the nonnegativity of [i follows. Now consider the sequence of sets 
'■— {C(2;) < Cmax ~ ?7n} where Tin is givcu by (|2.6[) . In view of (|2.7p . we have 

A+(V„) = ^-0, 

which clearly implies that supp/.t C E. 

By the compact embedding {Clip))* ^ (C°'^(2?))*, we deduce that fin ^ /j. strongly 
in (Co'^(X'))*. On the other hand, ([231) yields 



and the conclusion follows. 



Step 2. We now consider the general case. In view of the strong continuity of the functional 
under strong i?^- convergence, we can find a sequence {wn}neN C H^(T>; C) n C^(I?) 
such that for every n, 

\\Un- UtiWh^V) < £n, (2-17) 

and 

F,^{un) < F,^{un) + I , (2.18) 

so that sup„Fe^({t„) < +00. 

Given an arbitrary (p e Cq'^(I?) satisfying iVcpl < 1, we estimate 



{j{Un) ~ j{Un)) ■ V^ipdx 



V 

< \\Un - '«n||L2(X))||Vu„||i2(x,) + 1 1 || (-p) || V (-U„ - Un)\\L^(V) < CSnl ln£„| , 



using (|2.17p and (|2.18p together with Lemma 12.11 As a consequence, setting := 
a;~^curlj(u„), we have 

1 1 1^" -1^4 (cgi (15))- nZ^^ . (2.19) 

Applying Step 1 to {w„}, up to a subsequence we have /i„ — > ^ in (Cq'"^(I?))* for some 
nonnegative Radon measure fi supported by S. Then p.l9p yields — > in (Cq'^(I?))* 
and the proof is complete. □ 
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3. The lower bound inequality 

This section is devoted to the proof of claim (ii) in Theorem 11.11 that we summarize in 
the following result. 

Theorem 3.1. Assume that (HI) and (jl.2p hold. Let En — > 0"*" he an arbitrary sequence 
and let {unjneN C {!);£■) be such that 

fin — — curlj(M„) — > II strongly in (Cq'^(I?))* , (3-1) 

for some nonnegative Radon measure ii supported by S. Then, 

lim inf Fe„ (Un) > - Cmax ^(2?) • (3-2) 

n — *+oo 

In particular, if the left hand side in i3. is finite, then fi € H^^{T>). 

Proof. We will use in this proof the notations of the previous section. Without loss of 
generality, we may assume that 

liminf (m„) = lim F^^{un) < +oo . (3-3) 

n — »+oo n — »+C30 

Moreover, by Step 2 in the proof of Theorem 12. li we may also assume that {M„}„eN C 
H^{V; C) n C^iV). We shaU distinguish two cases. 

Case 1. We first assume that ln|lne„| < o(ti;„). We consider the family of vortex balls 
{-B"}ig7„ constructed in the proof of Proposition l2.2[ and we refer to it for the notations. 
Arguing as in (|2.3p we obtain 

> -7rAjlne„| -wJ;(//„,C) +o(l) (3.4) 



as n — > +00, where (•, •) denotes the duality pairing {Cq'^)*-Cq'^ . 
Combining the lower bound (|2.ip with p.4p . we infer that 

u;lF,^iun)>TTDn (ln(^^ ) - Ci] + ^ f |Vu„pda; + r!„ / W^C'jMdx 

> ttD,, (^In (^^^ - ^i) + ^ J^^ ^„ \\7un\^dx - iulifin, + o(l) , 

Since r„ = | ln£„|~'* and _D„ < 0{ujn) by Proposition l2.21 dividing the previous inequality 
by ujI yields 



FeSun)>7^l |Vu„pdx- (Ai„,C) +o(l) 

^ ^ / |Vu„pdx - CmaxM(25) + 0(1) ■ 

JV\VJBT 



In the last inequality, we have used (|3.ip and the fact that /i is supported by S. In view 
of claim (ii) in Proposition [^Hl we estimate 

/ \Wun\^dx> \u \j{un)\^dx. 

Jv\UB^ 1 + En JV,„\UB^ 
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Next we define 



so that by 



otherwise , 



0(1) > F.JUn) >\j^ \jn{x)\''dx - Cmax/i(2?) + o(l) . (3.6) 

Hence there exist a subsequence £„ ^ (not relabelled) and j* e i^(I?;M^) such that 
jn ^ j* weakly in L'^{'D; R^) as n — * +00. By lower semicontinuity, we have 

lim F,^{Ur,)>)-( |j;|2dT-Cmax/^(2?)- (3.7) 

It remains to tie the limit to the limit /i of the normalized weak Jacobians. To this aim 
we fix G !^{T>). Using Lemma [2.11 claim (i) in Proposition [5]T] and Holder Inequality, 
we estimate 



V ip ■ j{un) dx 
US" 



< (/:2(u,6/„Sr)) ' \\VipU\\Un\\LH'D)\\VUn\\L-(^V) 

<Cry2|ine„| =o(l). 
Since supp ip C for n large enough, we deduce that 

/ \7'^ip-j^dx— lim / \7-^ip-jndx= lim cj^^ / \7^{p ■ j{un) dx = 

= lim / W-'-ip ■ j{un) dx ^ - lim {pLn,(p)=- / (^sd/x. 

n — >+oo n — >+oo 

Consequently, 

curl = n in ^'(P) . (3.8) 

In particular, fi e H^^CD) since e L2(X';M2)^ 

Next we introduce E Hq{'D) to be the unique solution of 

\hf, = on dV . 



In view of p.Sp and the definition of h^, we have 

curl(> + V-^hf,) = in i?-i(r») , 
so that we can find € H^{'D) satisfying V.g^ = + V^/i^. Therefore, 







i 


Jv 


Jv 





since an integration by parts yields V^/i^ • Vg^ = (using the fact /i^ is constant on 
dV). Going back to p.7p . we infer that 



lim F,„(u„)>i/ iV/i^l^dx - Cmax/i(2?) • 



IV 

On the other hand, using the Green representation of h^^, we have 

'^h^\^ dx = 1- I I G{x,y)dn{x)dfj.{y) = . 
V ^ J Jvxv 
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and the conchision follows. 

Case 2. We now treat the case w„ < 0(ln|ln£„|). Consider the family of vortex balls 
{_B"}jgj^ (of size Sn = e~^^) constructed in Step 3 in the proof of Proposition 12.21 
Recall that this family satisfies Ds„ — Dn for n large. Combining the lower bound (|2.ip 
for the family {13"}j^j^ with (|3.4p . we derive 

^lF,M > (\n (^-^^ - ^i) + ^ l^^nfdx - iolif^nX) + 0(1) • 



Then arguing as in p.6p . we infer that 

1 
2 



Fe^{Un) > i / \jn{x)\^dx - CmaxM(25) + o(l) , 
V 



where 

'ij-^jiunix)) if X e Pe„ \ U,gJ„B7 
otherwise . 



As previously, up to a subsequence we have j„ ^ j* weakly in L^(I?;]R^), and 

lim Fs^{Un) > \ j \j*{x)\'^dx - CmaxM(2') • 
n-^+oo Z 

Now it remains to show that curlj* = /i in &'{T>), and then the proof can be completed 
as in Step 1. 

We proceed as before, taking an arbitrary Lp e 2>{'D) and using the weak formulation 
of the Jacobians. The key observation is that the contribution of the vortex balls will be 
negligible provided we can restrict our choice of test functions p> to functions constant 
in each vortex ball. This can be achieved thanks to |16[ Proposition 9.6], i.e., given an 
arbitrary ip g S){'D), there exists a modified function (pn which is constant on each ball 
and such that 

Il<y5 - '^n||co.°(X)) < C's.n^" , ||V(y5 - V(^„||l1(X)) < Cs„ (3.9) 

for each Q < a < \. Moreover, (pn has compact support in T) for n large enough. From 
and the L^-boundedness of the normalized currents jn , we derive 



/ (V^(p- V^^„) -jnda; < ||i„||i2(p)||V(^- V(^„|U2(P) < (3.10) 
Jv 

Then (|3.9p . p.lOp and the strong convergence of /i„ to ^ yield 

/ V'^p-j^dx— lim / V^p ■ jndx = lim / V^(pn ■ jndx = 

,-1 f V7±,~^ 



= lim W„ / V (pn-j{Un)dx^ lim - / Pnd^J.+ {^I- ^J.n,Pn) ] = (p dfl , 



and the conclusion follows. □ 



4. The upper bound inequality 

Throughout this section we shall use the following notation. For a nonnegative Radon 
measure ji E H^^{T>) compactly supported in I?, we denote by /i^ G Hq{T>) the solution 
of 

\h^=0 ondV. ^ ' ' 
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Using the Green representation h^{x) — G(x, y) d^{y), we have 
and the foUowing elementary lemma holds. 

Lemma 4.1. Let {fJ,n}neTi be a sequence of nonnegative Radon measure in H~^{D) with 
compact support in V. Assume that fin ^ fJ^ weakly* as measures on V as n ^ +00, for 
some n € H~^{T)) compactly supported in V. Then /(/^n) lilj) if and only ifh^^ 
strongly in H^{T>) as n +00. 

We may now start the proof of claim (Hi) in Theorem 11.11 in the case where the 
measure /i is absolutely continuous with respect to Ti'^LE. 

Proposition 4.1. Assume that (HI), (H2) and (|1.2p hold. Let Sn ^ 0+ be an arbitrary 
sequence. For every nonnegative Radon measure /i of the form 

fi = f{x)n'\-^ (4.2) 

with f e L°°(S), there exists a sequence {un}neN C H^('D;C) such that 

— curlj(u„) > /i strongly in {Cq'^{'D))* , 

UJn n—^+oo 



and 



lim (u„) = I{fJ,) - CmaxAi(2^) ■ 

n — >+oo 



Proof. Without loss of generality we may assume that / ^ 0, the case fi — being 
easily true. It is well known that a measure of the form (14. 2p belongs to H~^{'D), see e.g. 
[H Theorem 4.7.5]. 

We recall that the Dirichlet Green's function G in I? defined by (II. 3|) satisfies 



(i) G{x, y) > for every x E'D \ {y} and for every y G V; 
(a) for any compact set K CCT) there exists a constant Ck such that 



G{x,y) + ^ln|a:-y| 



< Ck (4.3) 



for all y E K and x E V. 



In the three first steps below, we assume that the density function / does not vanish 
on S, i.e., f > S Ti^-a.e. on E for some constant i5 > 0. The general case is considered in 
Step 4. 

Step 1. We will construct a trial function using the Green function G{x,y) in the spirit 
of [H]. Let 7 : [0,£] ^ E be an arclength parametrization of the curve E (so that 
£ = n\i:)), and define for t e [0,^], 

M{t) -.^ fi{j{[0,t])) . 

Then M(-) is strictly increasing, M{0) = 0, and M{e) = ^(E) = ^(P). Moreover M(-) is 
continuous since /i G H^^{'D) and thus atomless. 
Next we introduce for n large enough, 

_ [ cu„/i(E) " 
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where [•] denotes the integer part. Since M is continuous and increasing, we can define 
for fc = 0, . . .,!?„, 

Now we set for A: = 0, ... , D„, 



and we claim that 



Cu-^ <\al-a-,_,\<27rS-^u;-^ (4.4) 



for each k E {1, . . . , Z3„} and for some constant C > independent of n. Write Sfc := 
7([^fc-i,n, tk,n\) for fc = 1, . . . , Dn SO that Sfc is a smooth curve whose end-points are 
and a^k-i- construction, we have 

2w-i = Af (ife,„) - M{tk-i,n) = /^(Sfc) > Sn'{^k) > - Ofc-il . 

Now the curve S being smooth and {a^ — a^-il small for n large, we deduce 

Define the family of measures 

k—i 

We claim that fin ^ 1^ weakly* as measures on V and l^{T^)- To prove the 

weak* convergence of /2„, we fix an arbitrary function g Cq{T)). Observe that by (|4.4p 
and the smoothness of S, we have diam(Sfc) < Cquj~^ for a constant Cq independent of 
k and n. Therefore, using /i(Sfe) = 27rw~^ we derive 



where 



tpd^— ij; dfi. 



+ o(l) < C^i{V)osc{^,CoUJ-^) + o{l), 



osc(V',Cow„i) := sup |V'(a;) - V'(y)| 0, 

and the claim is proved. 

Next we must regularize the measure /2„ . Let us define for k — 1 , . . . , I?„ , 

fk,n{x) := jXS^„(a^) > 

where xb^ (a^) denotes the characteristic function of the ball B^^{a^). Set 

2 

/„(x) — ^/fc,„(a:) and fin '.^ fn{x)C^\_V . (4.5) 

^" fc=i 

Since e„ = o(a-',7^), the functions f^^n have disjoint supports for n large by (|4.4[) . As a 
consequence, fini'D) = fLniJ^) fJ-i^^)- Since /2„ /i, one may also easily check that 
fin ^ fJ- weakly* as measures on V. 

Step 2. According to (|4.ip . we introduce 
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Then 

/ \Vhn\^dx = Lol f I G(x, v) fn{x)fn{y) dxdy (4.6) 

= G(x,2/)/,,„(x)/,-„(y)dxdy. (4.7) 

We need to estimate the integral term in the right handside of (14. 6p . We proceed as in 
[T] and we provide some details for the reader's convenience. Let Ao CC I? be a small 
tubular neighborhood of S. Let < a < 1 be given small, and set Aq, := {{x,y) G 
V X V : \x — y\ < a}. Since G is continuous on (Ao x A/q) \ A^, and the support of //„ 
lies in Mq for n large, by the weak* convergence of jl„ to /i we have 



la ^ lim / / G{x, y) fn{x)fn{y)dxdy = 

i\A„ 

G{x, y)dii{x)dii{y) < ff G{x, y)dii{x)dii{y) = 21{p) . (4.8) 



Near the diagonal Aq, we split the sum in (|4.7p in two terms. Using (|4.3p we estimate 



D„ 



112 4^' E / / 2^) kn{x)kn{y)dxdy 



D„ 

^Xl// Gia'l +enZi,a'^ + enZ2)dzidz2 

/Si(0)xBi(0) 



1=1 

^4E/L„, ^,Jlln(-^)+c)dz,dz2 



1=1 



Si(0)xBi(0) 



<27rAjlne„|+0(w„), (4.9) 



and 



II 12 := 47r2 ^ / / G{x, y) /,,„(a:) fj,n{y)dxdy < 

0<\a^~a'^\<a " 

<C Yl H<~a]\\. (4.10) 

0<|a7-ay|<a 

By the smoothness of S and (|4.4p , there exists a constant cq > independent of n such 
that for every i j and every (x, y) € x Sj, 

|ar-a7l>co|x-y|. 

Since /i € II^^{T>) and it is supported by S, the map (x, y) ^ | ln(co|a:: — y\)\ belongs to 
L^(S X E, /i (g) /^). Therefore, by the Mean Value Theorem, for every i ^ j we can find a 
pair (a:", y") G x E^ such that 

|i // I ln(co|x - y\)\d^i{y)diJi{x) = | Hc^\x^ ~ y^\)\ , 
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noticing that 47r^cj„^ — fi{'E,i)^{Y^j). Applying the previous inequahty to (a;",?/") we 
deduce from (|4.10p that for n large, 

in:<c \HcoK-y^\)\ 

<CujI // \lii{co\x ~ y\)\d^J.{y)d^i{x) 

0<\a'l-aj\<a'' •'^i^^i 

<CloI(( \\n{c„\x ~ y\)\dti{y)dy.{x) . (4.11) 

■/■/(ExE)nA2Q 

On the other hand ^ E H^^{T>) so it is atomless, and thus fi (g) ^ does not charge 
{x = y} OV X T). Consequently the integral term in the right handside of (|4.1ip vanishes 
as a — > 0+. Hence, 

lim limsup uj:^'^in2 = . (4.12) 
Gathering (giH]), (gS]) and fil^ yields 

limsup I / / G{x, y) f„{x)fn{y) dxdy - 2'n:Dn\ ln£„| | < 

< limsup ( + limsup ((//^ - 2^i^„| ln£„|X2 + tj-^///") ) < 2/(^) , 

and in view of (|4.6p . we conclude that 
1 
2 



|V/i„pdx < TrDn I lne„| + Lull{fi) + o{iul) . (4.13) 



Step 3. We shall now define a complex- valued order parameter u„ associated to We 
proceed as follows. Since 

curl (-V-L/i„) = -A/i„ = w„/„ (4.14) 

is supported by Ufei3e„(a5J), we may locally define a phase 0„ in I? \ UfcSe„(a^) by 

V0„(a;) = -V-L/i„(x) for x e P \ UfcBe„ «) • 

In fact, since the balls i?£^(aj!) are pairwise disjoint (assuming n large enough) and the 
mass of oJnfn is quantized in each such ball, it is easy to show that is single-valued 
modulo 27r, i.e., for any smooth Jordan curve inside I? \ U/ci?e^ (aJJ), 



^ / V0„ • T e Z , 
/e 



27r 



where t : ^ is any smooth vector field tangent to 8. Hence exp(i0„(a;)) is well 
defined for every x £'D\ i)kBe„ i^k)- 

Then consider a smooth cut-off fimction p : K ^ [0, 1] such that p{t) = 1 for t > 2, 
and p(t) = for t < 1. Define 

{ ^ \ ^ ^ i^^n^ Jqj, g^j^g /c = 1, . . . , D„ , 

V £n y (4.15) 
1 otherwise , 

and observe that 

E,^ipr,,V) ^ 0{LUn) . (4.16) 
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Then set 

V„(a;)e'*"(") for x G 2? \ UfcB,„ (a^) , 
otherwise . 



Unix) := 



One may easily check that u„ e H^{T)\ C). We claim that 
A simple computation gives 



icurljK) ^ strongly in (4.17) 

n — *4-oo 



j(w„) = /O^V0„ = -plV^K a.e. in X* . 



Given ip e Cq' (I?) satisfying < 1, we deduce from (|4.14p and the previous indentity, 

1 



+ — / (a„ - 1)V(^ • V/i„dx . (4.18) 



JUfcB2.„(aJ) 

In view of the compact embedding {C^iV))* ^ (Cq'^(I?))*, /i„ — > ^ strongly in 
(Co'^(I?))*. Hence we can estimate using (|4.13p . 

lifin" f^,(p)\ < ||An - Ai||(c«-i(X))). +^n^('C^(UfcS2£„(a's!)))^^^||Vft.„||i2(x,) < 

< ||/i„ -/l||,„o,l.„^^. +Ce„|lne„|i/2 — , q, (4.19) 

and (|4.17p is proved. 

We now compute the energy F^^ (u„). We infer from ()4.13|) and (I4.16P that 



•D\UfcB,„(a^) 

<l [ |V/i„|' + 0{ujn) < nDn I lne„| + ujll{fi) + o{ljI) , (4.20) 
and it remains to evaluate the interaction with the rotation potential. First (|4.17p yields 



!^ / y^<:-j{un)dx-u,Mi^)u;i + o{cji). 

^Smax Jx> 



Arguing as in (|4.18p - (|4.19p we derive 

y-^C ■ jiun)dx = -w„(//„,C) + 0(e„w„|ln£„|^/2) 



V 

D. 



= -'^^Yl^ I C(a;)da; + o(£„|lne„p/2) = -27rL'„Cmax + o(e„|lne„|3/2), 

and consequently, 

I V^C- ]{un)dx^~^Dn\\nen\-U^.MV)ujl+o{Lol). (4.21) 
Combining ((i?^ with finally leads to 

Fe^ (U„) < I{p) - CmaxM(25) + o(l) • 
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In view of Theorem 13.11 the conclusion follows taking the limsup as n — > +00 in the 
previous inequality. 

Step 4- We now consider the case where the density / is allowed to vanish. Let {6k} C K. 
be a sequence decreasing to as fc ^ +00. Then for fc G N, we consider the measure 

Mfc :=A* + 4.H'LS= (/(a:)+4)HiLS. 

By monotone convergence, one has 

lifik) /(/i) and fik{V) ^liV). (4.22) 

k — *+oo k — >+oo 

Obviously also converges to ^ strongly in (Cq'^(I?))*. Applying Step 1 to Step 3, we 
find for every fc g N a sequence {u^}„gN C H^iV-.C) such that aj~^curl j(w^) — > Hk 
strongly in {C'^'^^iV))* and F^^{v^) Hpk) — CmaxA^fc(2') as n +00. Hence for every 
fc € N, we can find Nk G N such that for every n > Nk, 

IK'curlj(v^) - mII(c°-1(p))- < 2"'' 

and 

|i^,„(z;^)-/(A.) + C„axM2?)| <2-'=. 

Moreover we can assume without loss of generality that the sequence of integers {-/VfcjfcGN 
is strictly increasing. Therefore given any integer n large enough, there is a unique fc„ S N 
such that N^^ < n < N^^^i, and fc„ +00 as n ^ +00. We may then define m„ := w^" . 
By construction, the sequence {m„} satisfies the required properties. □ 

To consider the case of a general measure n in H^^{T>), we shall need the following 
continuity lemma. 

Lemma 4.2. Let fi be a nonnegative Radon measure in H~^(T>) such that supp fi CC T). 
For ^ G M^, let T^ji he the translated measure defined by 

T^niB) = fi{-^ + B) for any Borel set B C . 

Then there exists Q < 5 < dist (supp /i, dV) such that r^/i g H~^(T>) for every ^ G Bs{0), 
and the mapping ^ 1-^ h^:^ G H^{'D) is strongly continuous on Bs{0). 

Proof. For i5 > we set Vs := {x G , dist (a;, V) < 25}. Then choose 5>0 such that 
dT>s is smooth and 26 < dist (supp/z, dV). For every ^ G -85(0) we have V CC £, +T>s, 
supp T^fJ. — ^ + supp /i C I? and dist (i^ + supp /i, dV) > 6. 

Obviously /i G H^^{'Ds) and we can set h G H^{'Ds) to be the unique solution of 

{—/SJi = fi in Ds , 
h = on Vs. 

By our choice of 6, the function h is smooth in the ^-neighborhood of dV. Next, for 
^ G Bs{0) we denote by G H^{V) the function defined by := h{x — ^) for x G V. 
Observe that h^ G H^{V) and -Ah^ = r^fj, in V. Hence r^^ G H-^{V). 

Now consider a sequence {^„} C -65(0) such that ^„ ^ ^ G Bs{Q) as n — > +00. Denote 
hn := hr^^fj^. We have A{h^^—hn) = in I? and (h^^ — hn) = h{x—^n) on 91?. By standard 
elliptic estimates, (/i^^ — hn) strongly converges in H^{T>) to the harmonic function in 
V equal to h{x — ^) on 91?, that is h^ — /it^^. On the other hand, h^^ h^ strongly 
in H^(T>) by strong continuity of translations in H^. Therefore /i„ h^.^ strongly in 
H^{'D), and the proof is complete. □ 



22 S. Alama, L. Bronsard & V. Millot 



Theorem 4.1. Assume that (HI), (H2) and (|1.2p hold. Let e„ 0+ be an arbitrary 
sequence. For every nonnegative Radon measure ji G H^^(T>) supported by S, there exists 
a sequence {unjnGN C H^{'D;M.'^) such that 



and 



— curlj(u„) > /i strongly in {Cq'^ {T)))* , 



lim (u„) = /(^i) - CmaxAi(2') ■ 

n — ^+oo 



Proof. We shall prove that for any nonnegative Radon measure fi E H^^{T>) supported 
by E, there exists a sequence of nonnegative Radon measures {/^fcjfcgN of the form (|4.2p 
such that ^fe — > ^ strongly in (Cq'^(I?))*, /ife(2?) — > ^^(2^) and 

/(^fe) "> /(/i) asA:-^+oo. (4.23) 

Assuming that such a sequence exists, Proposition 14.11 yields for each k a sequence 
{v^}nm C iJi(X>;R2) such that uj-^cnrlj{v^) ^ifc strongly in {Cq^{V))* and 
-fe,i(^n) ^ ^(Mfe) ~ CmaxMfe(I') as 71 ^ +00. It then suffices to apply the diagonal argu- 
ment used in the proof of Proposition 14. 1[ Step 4, to construct the required sequence. 



Step 1. We first consider the case where S is a segment in V. Without loss of generality 
we may assume that E — [a, h] x {0} CC T) for some a, 6 G M with a < b. Assume in 
addition that E' :— supp^ CC]a, 6[x{0}. We shall regularize the measure using the 
following standard procedure. Consider a smooth function g G C°°(R) such that g > 0, 
supp £1 C [—1,1] and J^^g — 1. For a positive integer k and x = {xi,X2) S M^, we 
introduce gk{x) :— kg{kxi), and we define 

9k{x) j gk{x - y) dfi{y) . 

By construction, the function gk is nonnegative, smooth and supported by [a, b] x M. for 
k large enough. Next we define for k large the measure 

Mfc gkix)'H^\-T, . 

One may easily check that /i/j G H^^{T)), ^k{T^) m(^) and that pt^ ^ /i weakly* in 
the sense of measures on I? as fc ^ +00. In particular, /^^ — > strongly in [Cq'^{T>))* . 
We claim that (|4.23|) holds. Indeed, using Fubini's theorem we first derive that 

I{ti-k)^\ll G{z,z')dfik{z)dfik{z') = 



2 

1 
2 



(EnBi (x))x(SnSi (y)) 

Tc It 



G{z, z')gk{z~x)gk{z' -y) dHldHl, ) dfi{x)dfi{y) 



Next we observe that for k large enough, we have E n Bi{x) = {xi, 0) + Jk with Jk := 

(^, ■^) X {0} for every x = {xi, X2) G E'. Changing variables in (z, z') and using Fubini's 
theorem again, we obtain 

I{l^k) = \ll ill G{x + ty + e)gk{09k{e)dnldnl)d^i{x)dfi{y) 

^jjs'xs' yjjjkxjk / 

G{x + ^,y + C) dfi{x)dfi{y)^ QkiOBkie) dn\dn\, 
G{x, y) d{T^^i){x)d{T^.^Ji){y)] gk{0Qk{O dH\dn\, . 
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From the Green representation of /ire/^ we infer that for every Jk ^ Jk, 



G{x,y) d{T^n){x)d{T^, ^i){y) = / (V/i^^^) • (V/i^ dx . 
VxV Jv 

Then from Lemma 14.21 we deduce that the function 

0:(C,e')^ / i^K,^) ■ {Vhr^,^)dx 
Jv 

is continuous on -65(0) x Bs{Q) for some < S < dist (S, 91?). Therefore, 

hin /(/ifc) = lim Iff &{^,agk{0Qk{OdHldHl, = ^6(0,0) = , 
and (|4.23p is proved. 

Step 2. We now consider the case where S is a smooth embedded arc. We further assume 
that there exists a C^-diffeomorphism $ : I? ^ I? such that $(a;) = a; in a neighborhood 
of dD and E := is a segment compactly included in T). Let /i be a nonnegative 

Radon measure in H~^{'D) whose support is compactly included in the relative interior 
of S. Denote by p, the push-forward of fi through <I>, i.e., p, := ^#^J■. Then supp/2 is 
compactly included in the relative interior of S and /2 e H^^{V). Indeed, we easily check 
that 

m = \ ff_ v) dji{x)dm ^ \ ff G($(x), $(y)) d^l{x)d^l{y) 



< 



^ ff In - $(y)| dfi{x)dfi{y) + C 



< — / / In |x — y| dij{x)dii{y) + C < +00 . 



47r 



where we have used (|4.3p and the constant C only depends on S, /^(S) and || V$||loo(p). 

Therefore we can apply Step 1 to /2 to find a sequence of measures {p,k}keti of the 
form (|4.2[) such that supp/ife C E, /Ufc(2?) — + Pi'D), pk P- weakly* as measures on T), 
and /(yUfe) ^ /(/i) as fc ^ +00. Then we set fik ■— for every integer k. Observe 

that /Ltfc is of the form (|4.2p . Indeed, writing pk = 5fe(a;) LS with S C°(S), the area 
formula (see e.g. [S]) yields 

fik^gk0^x)\VrHx)\n^L^, 

where Vr$ denotes the tangential gradient of along E. Then one may check that 
IJ'ki'D) /i(2?), /ifc — ^ /i weakly* as measures on V as k ^ +00. 
We claim that (|4?251) holds. First write 

f f G{^-\x),<f-'{y)) dpk{x)dpk{y) 

'SxE 



S{<P~\x),<f-\y))dPkix)dpk{y) 

SxS 

+ T~ // In ( w X ^ w M ) dpk{x)dpk{y) 

= : 4- + //fe , (4.24) 
where S denotes the regular part of the Green function G, i.e., 

S{x, y) := G{x, y) + ^ In |a; - y\ 
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(which is a locaUy smooth function on V x V). Since fik converges weakly* as measures 
to /i, we have 

Mfc ® Mfe ^ M (8) A weakly* as measures on V x V , (4-25) 
and we deduce that 

\im h = l [[ S{^"\x),^-\y))dfl{x)djx{y) = \ II S{x,y) d^x{x)d^x{y) . 

(4.26) 

Next we consider a decreasing sequence q;„ — > 0. For every integer n, we introduce 
a smooth cut-off function Xn & C°°{'D x V) such that < x„ < 1, Xn{x,y) = if 
\^-^{x) - $~^(y)| > a„, and Xnix,y) = 1 if |$"i(a;) - $"^(y)| < a„+i. Note since 
fi G H~^{'D), the measure jl has no atoms, and hence /I® /I does not charge the diagonal 
{x = y} DD X v. Consequently, Xn ^ p,-a..e. vaV xV. Then write 

(4.27) 

By the choice of Xn and (I4.25p . we have for every n, 

.I'fJ''^^ = ij LJ' - ^"(-'2^))!'^ ( i^-i(.)-^-i(,)i ) • 

Next observe that 

< 1^ / ^ ,1 < forevery (x,y)el?xl?, (4.28) 

F-yl \^-^{x) - ^-^{y)\ \x-y\ 

for some constants Ci > and C2 > independent of a; and y. Since /(/i) < +00, 
estimate (|4.3I) tells us that the function In \x — y\ belongs to L^{T> xV^p,® p,). Therefore 
we may apply the dominated convergence theorem to derive 

Let us now treat the term Illji- We first deduce from (I4.28P that 
i /X _ Xn(a^,2/)ln (^^) dfik{x)dfik{y) < IIIJ: < 

<^ II Xnix,y)ln(-^)dflkix)dflk{y)- (4.30) 

Since the function 

G[x,y) 7y In I 1 =S{x,y) + In x - y In C^) 

ZTT \\x — y\J 2n ztt 

is locally smooth inV x V and I{pk) we infer from (|4.25p that 

, ^™ T~ // Xnix,y)ln( ] dflk{x)dflk{y) ^ 

^ Xnix,y)\n( ) dfl{x)dfi{y) fori = 1,2. (4.31) 
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Using that Xn — + /i (g) /i-a.e. and In \x — y\ belongs to L^{T> x T), ^ ® ji), we infer as 
previously that 

lim Xn{x,y)\n{-^\dii{x)dn{y)^Q for z = 1,2. (4.32) 

Combining (li^O)) . (|OT1) and (li?^ we derive 

lim lim inf 1 11^ = lim lim sup III^ = , 

which yields together with (|4.27p and (|4.29p . 

Then follows gathering g^, (|4:26l) and g^S]). 

S'tep 5. We now consider the general S case. If E is an embedded arc, we may assume 
without loss of generality that E C E' for some C^-Jordan curve E' compactly included 
in T). Hence it suffices to consider the case where E is a Jordan curve. We shall use the 
following lemma. Its proof is postponed at the end of the section. 

Lemma 4.3. Assume that Yi is a C^-Jordan curve. Then there exists Si > such that 
for every xq G E, there exists a -diffeomorphism ^ : D —> D satisfying $(x) = x in a 
neighborhood of dT> and such that <I>(En Bs-^{xq)) is a segment compactly included in T). 

Now let 7 : [0, 1] ^ E be a constant speed paramctrization of E. Let be a positive 
integer to be chosen and set tn = n/N for n = 0, . . . , N, and 

E„ := 7([^n-i,i«]) for n = 1, . . . ,7V. 

We choose N in such a way that diam(E„) < Si for each n, where the constant Si is given 
by Lemma 14.31 Setting x„ — 7((tn-i + tn)/2) for 7i = 1, . . . , A'^, we can apply Lemma 14.31 
to each Xn to find a C^-diffeomorphism $„ : I? ^ I? such that <I>„(E„) is a segment 
compactly included in V, and ^^(x) = a; in a neighborhood of dT>. 

Let /i be an arbitrary nonnegative Radon measure in H^^{'D) supported by E. Con- 
sider a decreasing sequence — s- and define for k large enough, 

Yn-^li[tn-i + ak,tn-ak]), AiJi:=A*LEf; for n 1, . . . , iV. 

Oviously e H^^CD) with supp^^ C E^. Applying Step 2 for each n and k, we 
find a sequence of measures {MnmjmeN of the form (14. 2|) such that supp/ij^ „j C E„, 
fJ-tmC^) ^^(2?), lJ-n,m ^ l^n wcakly* as measures on V, and I{p'^^^) as 
m +00. Define the measures 

N N 

n— 1 n—1 

Then ^^iV) — > fJ-^iJ^) and /i^^ ^ /i*^ weakly* as measures on I? as m ^ +00. In addition, 
from Lemma l4.1l we infer that h^h h^k strongly in H^{V) for every integers n and 
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k. Hence 

N „ „ ^ N 



G{x,y)d^i'ld^,';^I{^/) 

V 



Next recall that /i is atomless. Hence, by monotone convergence we have ^ki'D) 
IJl{T>) and /(/^'^) as fc ^ +00, as well as the weak* convergence of to /x. 

Consequently, 

hm hni \^ll{V)-^x{V)\= Win lim ~ 

= lim lim ll/i^, - mII(cO-1(X)))* = 

(here we use again the compact embedding (Cq(I?))* ^ (Cg'^(I?))*), and the conclusion 
follows for a suitable diagonal sequence = f^m ■ n 



Proof of Lemma \4.S\ By assumption on E, there exists 60 > such that for every xq & S, 
S n -6250(2^0) is the graph of a -function and -B2i5o(^o) <^ ^- Now fix G S and write 
every x G D a.s x = Xq + st + ir"*" where r denotes a unit tangent vector to S at a;o- 

Then E n B2So{xo) = {xq + ST + /(s)r-'- , S e [Smin, Smax]} for some > Snun > -2(5o, 

< Smax < 26q, and a C^-function / : [smin, Smax] ^ R satisfying /(O) = /'(O) = 0. Since 
E is C^, there exists a constant k > which only depends on S such that |/"(s)| < k 
for every s € [Smin, Smax]- Hence we may choose Sq smaller if necessary (uniformly with 
respect to xq) in such a way that |/'| < 1. Then Smin < —Sq, Smax > Sq and E n Bs{xq) 
is still a connected arc for any d < 2So. 

Set Si :— (5o/(2 + k). We claim that E n Bsi (xq) satisfies the requirement. Indeed, we 
may construct a C^-diffeomorphism $ : P — > I? as follows. Consider a smooth cut-off 
function x : ^ M such that < x < 1, x{x) = 1 if \x\ < Si, x{x) = if \x\ > 2Sq and 
|Vx| < <5(7^ Then we set for x eV, 



^{x) X - x{x - xo)f{{x - Xq) 



T)T 



± 



The reader may check that $ maps V into T), $ is one-to-one and defines a C^- 
diffeomorphism. Moreover $(E n i?5^(a;o)) = {xq + st , —Si < s < Si} is a segment 
compactly included in V. □ 

Proof of Corollary 11.11 Step 1. For any nonnegative Radon measure /i supported by 
E we have 

/(m) - U.M'D) > h {^{V)f - CmaxM(25) , (4.34) 

and equality holds if and only ii ^ = Xji^, for some constant A > 0. We recall that /x* 
is the unique minimizer of / among all probability measures supported by E and that 
/* :— I{fJ,^,). The existence and uniqueness of fl^, is classical, and we refer to [T3] for 
further details. Optimizing (|4.34p with respect to A for measures of the form /i — A/x* , we 
derive that ^^/x* is the unique minimizer of /x 1-^ ~ Cnax/^(2?) over all nonnegative 
Radon measures supported by E. 

Step 2. Let e„ 0"*" be an arbitrary sequence. The existence of a minimizer ?x„ of 
is classical and follows from standard arguments based on coercivity and lower semi- 
continuity properties of F^^. We first observe that Fe^{un) < ^6„(1) = 0. Hence, by 
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Theorem II. 1[ there exists a subsequence {e„j. } such that 

j("«J f^o strongly in (Cq'^CD))* as fc — > +00, 

for some nonnegative Radon measure /ip G H^^(T>) supported by S. Moreover, 

hminf ^^e„, ("«J > /(a^o) - CmaxMo(2?) • (4.35) 

On the other hand, by Theorem II. li any nonnegative Radon measure ^ S H~^(T>) 
supported by E can be strongly approximated in (Cq'^(I?))* by some sequence {ij^nlK^k)} 
with {vk] C H^{T>\C) satisfying 

lim Fe„ (Wfe) = /(^) - CmaxM(2?) ■ 

— ^+00 ^ 

Since Fg^^{u„^) < F^^^{vk) we infer that /ig minimizes /i 1— > /(/i) — Qmax^-i'D) over all 
nonnegative Radon measures supported by S. Consequently, /xq = ^ff^fi* and the liminf 
in (|4.35p is actually a limit. Then the result along the full sequence {£«} follows from a 
standard argument on the uniqueness of the limit. □ 

5. r convergence analysis for annular domains 

In this section we briefly show how to extend the above techniques to the case of a 
multiply connected domain. The method we outline here may be applied for any finite 
number of holes (see [2]), but for simplicity we restrict to domains which are topological 
annuli. Let V denote a simply connected domain in with smooth boundary, and 
B CC V a smooth, simply connected domain compactly contained inside 2?. Then let 
A:^V\B. For u e H^{A; C) we define the functional 




Here the given vector field V : ^ is assumed (for simplicity) to be locally Lipschitz 
continuous. We are interested in the asymptotic behavior of as e — > 0, with an angular 
speed as in (|1.2p . 

5.1. Asymptotic vorticity of the hole 

For multiply connected domains, the highest order term in an expansion of the minimal 
energy is partially due to the turning of the phase of a minimizer around the holes. The 
first step in studying vortices in the interior is to identify the asymptotic vorticity of the 
hole, and then split the energy into contributions from the hole and from the interior. To 
this purpose we first study the minimization of the functional over S^'^-valued maps. 
Observe that for §^-valued maps, the functional only depends on the angular speed 
and not anymore on e itself, i.e., for every u G H^{A;§'^), 

Uu)^nn(u) -.^ j^\^\^u\''~nV{x)-j{u)^ dx, 

with $7 = rig. We are therefore interested in minimizing Tin over the class iJ^(^;§^), 
and here > could be any positive parameter. It is well known that maps in H^{A; S^) 
are classified by their topological degree, i.e., their winding number around the hole B. 
Hence, minimizing first in each homotopy class and then choosing the lowest energy level, 
one reaches the minimum of the energy of Tin, i.e., 

minTio = imiig{d, il) , (5-1) 
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where 

g{d,n) ■.= ijim{nn{u) : u e H\A;S^) , degu = d}. (5.2) 
Concerning the minimization problem (|5.2p . we have the foUowing result. 

Proposition 5.1. For every d G Z. the minimization problem (j5.2p admits a unique 
solution up to a (complex) multiplicative constant of modulus one. Moreover, 

g{d, {iV'fdP - ^^'IV^I'} dx , (5.3) 

where $d is the unique solution of the linear equation 

'-A$rf = ricurlF in A, 

$d = ondV, 

$d = const. on dB , (5.4) 

^^'^ 2Trd-n I V-T 



dB JOB 

Proof. We follow here some of the arguments in |6, Chap. 1], and we provide some 
details for the reader conveniance. 

Step 1. We claim that for any u E H^{A;§>^) such that degw = d, we have 
nn{u) >\ I {|V$<ip-f]2|y|2} 

^ J A 

Indeed, we first observe that curlj(u) = since u is §^-valued. On the other hand, 
V^$c( + VlV is also curl- free and 

JdB JdB C'J^ Jqb 

so that we can find a scalar function H G H^{A) such that j{u) = \7H + V^<I>d + flV. 
Since u is §^-valued, we have jVup = |j(u)p, and thus 

Hn{u) = ^lj\j{u)-nV\'-n'\V\^} dx 

= \ I {|V$<j|2-r!2|T/p} dx+]- f \VH\^dx+ f V^^d-VHdx. 
^ Ja ^ Ja J a 

Then, using the fact that the function $rf is constant on dA, an integration by parts 

yields jj^V^^d ■ VH dx = and the claim follows. 

Step 2. We claim that there exists Ud £ H^{A; such that degUd — d and 
Indeed, since curlCV^^d + VtV) = and 



r / (v^*d + nv)-T ^d£ 

JdB 
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we may locally define a scalar function ■0 in A such that 

VV' = v^^d + ^v . 

Then Ud '■= exp(zi/') is well defined and satisfies the required properties. Clearly the 
construction of Ud is unique modulo a constant phase, and the proof is complete. □ 
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In order to solve problem (|5.ip , it now suffices to express (j5.3p explicitely in terms of 
the integer d. To this purpose, we first introduce the solution ^ of the linear problem 

= in A, 

^ = on 9X» , (5.5) 
_^=1 ondB. 

The function ^ is smooth in A and < ^ < 1 by the maximum principle. Moreover, 
the Dirichlet energy of ^ is the so-called i7^-capacity of B inside V which we denote by 
cap(S), i.e., 

cap(S) := / \Vi\'^dx = ~ f 1^ >0. (5.6) 

Next we consider the unique solution C, of 

f-AC = curly in^, 
I C = on oA . 

From the Lipschitz assumption on V and standard elliptic regularity, we infer that C 
belongs to Cq'°'{A) for every < a < 1. We set 



dV 



d 



V 



Observing that (15. 7p implies 



dB JdB 

we find that for every integer d, the function determined by (|5.4p is explicitly given 
by 

(5.8) 

V cap(S) ^ 
Moreover, using (|5.5p . (|5.6p and (|5.7p we readily obtain that for every d e Z, 

As a consequence, an integer dn is a minimizer in (|5.ip if and only if do minimizes the 
function d G Z Itv^^ — 27r(i|. 

We may now state our result concerning problem (|5.ip . 

Theorem 5.1. Up to multiplicative constants of modulus one, the minimization problem 
(|5.ip admits exactly two solutions (of distinct topological degree) if jv^/t^ is an odd 
integer, and a unique solution otherwise. Moreover, if da & is a minimizer in (|5.ip . 
then da G {[-^57^]; + 1} where [■] denotes the integer part, and 

minHn = / ~ |VCP} da; + C»(l) as n ^ +00 . (5.9) 

^ J A 
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5.2. The T convergence result 

To state the parallel F-convergence result for the anular domain case we must give more 
specific hypotheses on the potential V and the angular speed fl^. In addition to the 
Lipschitz regularity, we assume in the sequel that V satisfies the following assumptions: 

(HI') the solution ( of (|5.7p is such that (^max ■— niax^^^ ({x) — max^^^ \C{x)\ > 0; 

(H2') the set I] :— {x ^ A : C{x) — Cmax} CC ^ is a Jordan curve or a simple embedded 
arc of class C^. 

We note that in (Hl')j the assumption that Cmax is achieved at positive values of ( is not 
restrictive. Indeed, considering the complex conjugate of an admissible function replaces 
y by — y in the energy and hence C by — C. 

As for the simply connected domain case, we assume that fi^ is near the critical value 
needed for the presence of vortices which again reads 

n,^}^+u;{s), (5.10) 

for some positive function lu satisfying uj{s) — > +oo with a'(e) < o(|lne|) as e — > 0^, 
exactly as in l|1.2p . 

In the sequel, for an arbitrary sequence ^ 0"*", we will denote by u* a minimizer 
of Wn„, «-e., a solution of (|5.ip . and its corresponding topological degree will be denoted 
by dn- For brievety we shall also write (|5.8p as 

*n := = QnC + ^nC With «„ := ^ — . 

cap(o) 

We emphasize that = 0(1) as n — > +oo thanks to Theorem 15. II 
For V £ H^{A; C), we now define 

The following proposition shows that the functional F^^^ {Un^) captures the energy in- 
duced by interior vorticity of a given configuration u. 

Proposition 5.2 (Energy decomposition). Assume that (|5.10p holds. Let e„ 0+ 
be an arbitrary sequence and {un}neri C H^(A; C). Then, setting w„ — u^Un G H^{A; C), 
we have 

sup ri~^Je^(u„) < +00 if and only if sup il~'^uj'^F^^{vn) < +oo ■ (5-11) 

n n 

Moreover, if one of the conditions in (j5.1ip holds, then 

Js^{un) ^ minHn,, + ^^nFe„{vn) + o{l) asn^+oo. (5.12) 

Proof. Straightforward computations yield 

^ = ^ + 7i +j{Un)-J[Vn), 



and 

j{un) = |u„pj(0 + JM ■ (5.13) 
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By the proof of Proposition (STTJ we have = V^$„ + VlnV . Hence, 



Then we observe that 
Therefore, 

Since — |w„|, we also have the same estimate as above with Eg^{vn^A) instead of 
Eg^{umA). Assuming that one of the conditions in (jS.lip holds and arguing as in 



we derive that either Ei.^^{un,A) < 0(ri„), or E^^ (t;„ , A) < 0{Qn)- Consequently, 
if one of the conditions in (|5.1ip is satisfied, (|5.12p holds and the conclusion follows 
combining (|5.12p with ([5J| . □ 

For a nonnegative Radon measure on A, we define 



■= \ 1 1 G[x, y) dfi{x) dfi{y) , 
^ ' JaxA 



where the function G denotes the Dirichlct Green's function of the domain A, i.e., for 
every y £ A, (?(•, y) is the solution of 

i-AG{;y)^Sy mS)'{A), 
|G(-,?/) ==0 ondA. 

We may now state the F-convergence result for annular domains which involves the family 
of "reduced" functionals {i^gj^^o- 

Theorem 5.2. Assume that (HI'), (H2') and (|5.10p hold. Let £„ 0+ be an arbitrary 
sequence. Then, 

(i) for any C H^{'D;'C) satisfying sup„_F£^(i;„) < +oo. there exist a subse- 

quence (not relabelled) and a nonnegative Radon measure ^ in H~^(A) supported by 
E such that 

curl j(w„) — > /i strongly in {Gq'^ (A))* ; (5.15) 



n — >+oo 



w(en) 

(ii) for any {vn}neN C H^(T);C) such that (j5.15p holds for some nonnegative Radon 
measure ^ in H^^(A) supported by E, we have 

liminf F^^ («„) > - Cmax Ai(.4.) ; 

n — >+C30 

(Hi) for any nonnegative Radon measure /i in H^^{A) supported by S, there exists a 
sequence {vn}neti C H^{A;C) such that (j5.15p holds and 

lim Fe„ (u„) = ~ Cmax K-^) ■ 

n — ^+oo 

As in the simply connected case this F-convergence result could lead to the asymptotic 
description of the vorticity in i^^-global minimizers. Actually Theorem [5T2] combined with 
Proposition 15.21 also gives the asymptotic behavior of vorticity in J^-global minimizers. 
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The key observation here is that niinimizers for yield quasi- minimizers for F,,, and 
conversely. 

Corollary 5.1. Assume that (HI'), (H2') and (|5.10p hold. Let En 0+ be an arbitrary 
sequence. For every integer n g N, let u„ € H^{A;'C) be a minimizer of J^^^{-). Then, 

1 , \ Cmax _ J , . /^O.l/ 



■curlj(u„) > -:t^M* strongly in {Cq (A))* 



where /I* is the unique minimizer of !{■) over all probability measures supported by S, 
and := /(/!*). In addition, 

JeM = -^Jj\y\' ' IVCP} dx - ^^col + o{ul) . (5.16) 
We conclude this subsection with an elementary example motivated by [11213] . 

Example 5.1. Assume that V = Bi{0), B = Bp{0) for some < p < R and V{x) = a;-"-. 
Then the solution C of ()5.7p is given by 

, ~ , , , p'^lnR- R^lnp 

^(")='— +2h^^"'"'+ 2HR/P) ■ 

In particular, the set E is given by the concentric circle B^, (0) with 



^ e(p,i?). 



2\u{R/p) 



Here again, the uniform measure p^ — (27rr*)^^d?i^ LS turns out to be the Green 
equilibrium measure for E in A, i.e., I{p*) = /*. Indeed, one may easily check that the 
function 

ln(i?/n) 



/i*(a;) 



27r(ln(i?/r,) + ln(n/p)) 
ln(n/p) 



ln(|x|//9) if p < \x\ < , 
ln(i?/|x|) ifr,<\x\<R, 



27r(ln(i?/r,) + ln(r*/p)) 

solves —Ah.^, = /i* in A with = 0. Hence h^,{x) — jj^G{x,y) dpi,{y), and since /i, 

is constant on S the conclusion follows from Theorem II. 5. 12 in [T4]. 



5.3. Compactness of normalised weak Jacobians 

In this subsection we shall be concerned with the proof of claim (i) in Theorem 15.21 
We consider an arbitrary sequence £„ 0+. For any {w„}„gN C H^iV;C) satisfying 
sup„ Fs„{vn) < +00, we first derive exactly as in Lemma [XT] the estimates 

^e„(««,^) < 0(|lne„n and \\vjLHA)<Oil). (5.17) 

Hence, assuming in addition that {wnlneN C C^{A), we can apply the vortex ball con- 
struction in Proposition 12.1! with A in place of V, and Ae {x E A : dist(x, dA) > e} 
in place of T)^. We choose again r — rn := |lne„|~^, thus obtaining a finite collection 
of disjoint closed balls {i?(a", Pi,n)}ie/„ (written :— B{a2 , pi^n)), with associated 
degrees di^n and total approximate vorticity 

as in Proposition [2T1J 
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The first difference with the simply connected case arises in an estimate analogue to 
(|2.3p , as we must take into account an additional contribution to the potential term due 
to the boundary dB. Since one or more of the vortex balls may intersect dAe \ dV^, 
we will need to perturb this boundary slightly. When calculating the boundary term it 
will be convenient to choose a level set of . For < s < i < 1, denote by 

at -.^ {x eA: £,{x) ^ t}, and As^t := {x e A : s < < t} . (5.18) 

As an easy consequence of the Maximum Principle and the Hopf boundary lemma, each 
curve (Tj is smooth and the family {crt}o<t<i realizes a foliation of A. Then, for every 
t e (0, 1) the curve at is diffeomorphic to dB, and the set At^i is a neighborhood of dB 
in A. Now we shall choose an appropriate level set of ^. Define 

Jn = {te (0, 1) : atn (U,e/„Bn = 0} . (5.19) 

We note that the measure of the complement (0, 1) \ J„ is of the same order as r„ = 
I ln£:„|~^. Hence we can find t,i G J„ such that the level curve 7„ := at^ satisfies 

e„ < dist (7„,aB) < 0(|lne„r3), 

and consequently i„ = 1 + 0{\ lne„|~^). 

This construction allows us to define to topological degree of v„ around 7„ since w„ 
does not vanish on 7„, i.e., 

6n ■■= deg( f^,7n ) • 

Then an approximate total vorticity in A of the configuration Vn is given by |<5„| + Dn- 
We may now state the following proposition which parallels Proposition 12.21 

Proposition 5.3. Assume that (HI') and ((5TTU)) hold. Let {vn}neN C H^iV^QnC^iA), 
Dn and (5„ as above. Then \6n \ + Dn < 0(a;„). 

Proof. Step 1. Arguing exactly as in the proof of (|2.3p . we first derive 

I V^C • j{vn) = -27rf]„ d^.nCiaf) + o(l) . (5.20) 

Then following essentially the proof of Lemma 3.4 in [1], we obtain 

/ V^^J(«n) = -2^ Vd,,„C(an-2^i„J„ + o(l). (5.21) 

Here the fact that 7„ is a level set of ^ is essential in obtaining the degree (5„ from the 
boundary term when integrating by parts. 

Ste-p 2. As in the proof of Proposition 12. 2i we may assume that LOn < 0{Dn + \5n\)- From 
(I5.17P and claim (iv) in Proposition 12.11 we have D„ < 0(|ln£„|). Then we infer from 
(1231), (jog)) and ([5:iT|) that 

O(iol) > LolFeAvn) > 7ri?„ (| In £„ | ~ Cln|lne„|) - 27rr!„ ^ d^^nCK) - 2Tr\an\Dn 

-2T:\an\\5n\+ [ iVvnf dx + 0(1) , (5.22) 

JA\Uiei„B^ 

where we have used the fact that < ^ < 1 and < t„ < 1. 
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Next we consider a sequence 7y„ — > as in (12. 6p , and we group the vortex balls into 
different classes as in the proof of Proposition !^?^ fwe refer to it for the notation). Exactly 
as in (pTTOl and ((2?T2)) . we derive that 

TrD*Jlnen\ - 2tt^I^ ^ dunCiaf) > -27rc^„Cmaxi?: , (5.23) 

and 

nD+\\iien\ - 27rr!„ ^ d,,r.CK) > Cn^rinD+ . (5.24) 

For negative degrees, we observe that {C < 0} n E = since ^max = ICImax > 0. Hence 
we can estimate as for the class /+, 

> Cf}„r?„Z?;^ . (5.25) 

Inserting ([g?^ . (j??^ and ([g?^ in (lO^ yields 

O{iol) > -nCiDn In | ln£„| - 27rcj„CmaxAi + C?7„f^„(i:»+ + D") 

-27r|a„||,5„| + i / \Vvn\''dx. (5.26) 

Using the fact that |a„| = 0(1), we easily deduce the estimate 

' ^'-''^ 

for a constant C > independent of n. 
We claim that 

/ \Vvn\^dx>C{Dn + \Sn\f . (5.28) 

JA\UB^ 

Accepting (|5.28p . we infer from (|5.26p that 

(Dn + \Sn\ f - Cmax{cc;„, In | ln£„|}(i?„ + |5„|) < ©(cjfj , 
which clearly implies 

Dn + \Sn\ < 0(max{cc;„,ln| ln£„|}) . (5.29) 

To prove (|5.28p we introduce 
to := 1/2 min{t G (0, 1) : cii n E 7^ 0} , ti := 1/2(1 + max{t € (0, 1) : CTj n E ^ 0}) , 
and 

{t € (e„,io) : <Jt n {i)^eI„B^) = 0} , {t G (ii,i„) : ct* n (U,e/„Sn = 0} . 

Then dist (crt, E) > C > for every t & J^U and min(£i(ji), £i(-^n)) > C > for a 
constant C independent of n. Next we consider 

Dnit) := deg ( -^,crt 



If 2\Sn\ > Dn and n large enough, we estimate using (I5.27p . 

\Dn{t)\>\\Sn\-D+-D-\>^\Sn\>C{Dn + \5n\) for cvery t G . (5.30) 
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In the opposite case 2|(5„| < £>„ (and n large), we have in view of (|5.27p . 

\Dr,{t)\>\Dn~D+-2D-^\6,,\\>^D.^>C{Dr, + \5n\) for every i e ^ . (5.31) 

Set Vn '■= Vn/\vn\- Using claim (ii) in Proposition 12.11 the Coarea Formula and Jensen 
Inequality, we derive from (|5.30p and (|5.3ip that 




>C j ( f \vr.AVrVn\'dn')dt>C f dt > C {D,, + \5.r,\f , 

and (|5.28p is proved (here we have also used the fact that |V^| does not vanish in A). 

Step 3. If In I In e„ I = o(ti;„), the conclusion follows from (|5.27p and (|5.29p . If cj„ = 
0(ln I lne„|), we refine the lower bound by growing the vortex balls as already performed 
in Step 3 of the proof of Proposition 12.21 Following the same arguments with minor 
modifications yields the announced result, so we omit the details. □ 

Proof of Theorem 15.21 claim (i). In view of Proposition 15.31 and estimate (|5.27p . we 
can follow the proof of Theorem 12. 11 considering first a sequence {wnjneN C H^(T>; C) fl 
C^{A) and then the general case. □ 



5.4. The lower bound inequality 

Proof of Theorem l5.2L claim.(ii). We shall use the notations of Subsection l5.3l With- 
out loss of generality, we may assume that 

liminf Fe„(t;„) = lim Fe„ (t;„) < +oo . (5.32) 

n — *+oo n — *+oo 

As in the proof of Theorem l3.11 we may also assume that {wnjnsN C H^{T>; C) D C^(2?). 
We shall only consider the case In | lne„| = o(w„) since the other case can be completed 
as we already pursued in the proof of Theorem 13.11 
By Proposition 15 . 31 and (|5.2ip . we have 



V^e-j(v«)da; =0(w„). 

A 

Then we can argue exactly as in p.6p to derive that 



FeAVr.) >^f \j{Vn)\^ ~ CmaxA^(^) + o(l) , 



Setting jn (x) as in (|3.5p (with Ae^ in place of ) , up to a subsequence we have j„ ^ j* 
weakly in L^(.4;R^) as n +oo. By lower semicontinuity, we obtain 

liminf Fe„(w„) > i / - Cmax/i(^) • (5.33) 

In addition, arguing as in the proof of (|3.8p . we deduce that 

curl = ^ in ^'{A) , 
and thus fi G H-'^(V) since G L'^{A;R'^). 
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Next we introduce /i^ e i/p {A) to be the unique solution of 



and we set 



/im = on dA , 



By construction, we have curl (j* + V^/i^) = in H^^{A) and Jg^ij* + ^'^h^) • r = 0. 

Hence there exists 5 € -^^(-4) such that j» + S/-^h^ = \7g. Arguing as in the proof of 
Proposition l5.ll we derive 

J*\'= I \^K?+ I |Vgp> / \yh,\\ (5.34) 

A J A J A J A 

Then using (|5.5p and h^^^j^ = 0, we obtain V/i^ • = 0, so that 

Finally, using the Green representation of /i^ we obtain 

l f iVh^fdx^l ff G{x,v)dii{x)dii{v), (5.36) 

and the conclusion follows gathering (|5.33p . (|5.34p . (|5.35p and (|5.36p . □ 



5.5. The upper bound inequality 

Proof of Theorem I5.2|, claim ('iiij. We present here the proof in the case where the 
measure fi G H^^{A) is absolutely continuous with respect to LS, and more precisely 
for /i of the form (j4.2p with a nonvanishing density function /. The general case follows 
by approximation as already pursued in Section ID 

For such a measure fj, we first proceed exactly as in Step 1 of the proof of Proposi- 
tion [¥TT1 and we refer to it for the notation. For each integer n, we consider the function 
/„ as defined in (|4.5p and we set /t„ := fnC^\—A. Then fin ^ n weakly* as measures on 
^ as n ^ +00. Next we consider the solution ft,„ of 

-Ahn^UJnfn in A , 

hn — on dA . 

Arguing exactly as in Step 2 of the proof of Proposition 14.11 we derive that 

i / |V/i„|2 dx < ^AJ lne„| + c^^/(m) + o{ujI) . 
^ J A 

Next we introduce 

j_ r dhn] 
271" JdB dv J 

where [•] denotes the integer part. Noticing that k„ = 0(1) we deduce 

<7r^„|ln£„|+w^J(M) + o(w^). (5.37) 



cap(;B) 



^ with Kr, 



OB dv 



27r 
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In view of (j5.6p we have 

and since is harmonic in A, 

1 f dhn 1 f dhn 



r dhn _ 


' }_ r dhn 


JdB 





= — 1 for every k = 1 , • ■ • , Dn 



Hence, for any smooth Jordan curve 8 inside ^\ U/ci3e„(a^), 

where r : O ^ §^ is any smooth vector field tangent to 0. Consequently, we may locally 
define a phase 0„ in ^ \ DkBe„ (a^ ) by 

V(t)n{x) = -V^hnix), 2;e^\UfcSe„(a^), 

and then the map exp(i0„(a;)) is well defined for every x £ A \ UfcBe„(a^). 
Finally we consider a profile function p„ as defined in (|4.15p and we set 



Vnix) 



Vn(2;)e'*"(") forxe^\UfeS,„(a«), 
otherwise . 



Since — A/i„ = LOnfj'n in A, using ()5.37|) we may proceed as in the proof of Proposition l4.II 
Step 3, to prove that 

w^^'^curl j(wn) — > /i strongly in {Cq'^{A))* as n ^ +oo , 

and that 

E,^{vn,A) < TTDn\\nen\ + ujllifi) + o{ujI) . (5.38) 
To evaluate the rotation part of the energy, we first argue as for (|4.2ip to obtain 

/ VC^ • j{vn) dx = -ttAiI ln£„| - U^MA)ujI + o{ujI) . (5.39) 

J A 

In a similar way we derive that 

an / V-^^ • j{vn) dx = -On / VK ■VCdx + o(l) = -anKn + o(l) = 0(1) , (5.40) 
J A J A 

since a„ — 0{l) and /t„ — 0{l). Then the conclusion follows gathering (|5.38p . (|5.39p 
and (jOO)) . □ 



5.6. Application to - minimizers 

Proof of Corollav lS.li Let £„ ^ 0"^ be an arbitrary sequence. As for CoroUarv ll.li the 

existence of minimizcrs for J^^ and F^^ is classical, we omit the details. Then for every 
rt e N, let Un E H^{A;C) be a minimizer of Je„. Since J^^{un) < minTisi,^ = O(ri^), we 
can apply Proposition 15 . 21 to the sequence {wri}„gN to infer that 

min7io„ > JeAun) = min7io„ + ujlF-^^Vn) + o(f ) , (5.41) 

with Vn ■= UnUn- On the other hand, for any minimizer u„ £ H^{A; C) of Fe„, we have 
Fs^i^n) < -F'e„(l) = 0. Hence Proposition 15.21 vields 

Js„ {un) < Je„ «w„) = minWo,, + ujI min + o(l) . (5.42) 
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Combining (|5.4ip with (|5.42p . we infer that oj^F^^ (u„) < o(l) and that is a sequence 
of quasi-minimizcrs for {i^g^}, i.e., 

FeAvn) = minF,„ + o(l) (5.43) 

as n ^ +00. By Theorem 15.21 there exists a subsequence {Snk} such that 

^n^j i^n,) -> Mo 

strongly in {Cq"^{A))* for some nonnegative Radon measure hq € H~^{A) supprted by 
S. Using (|5.43p together with claims (ii) and (Hi) in Theorem 15.21 we deduce as in the 
proof of CoroUarv l 1 . 1 1 that /io minimizes /i /(/i) — CmaxM(-^) over all nonnegative Radon 
measures supported by E. The minimizer is again unique and given by 

Cm ax _ 

whence the convergence of uj^^j{vn) along the full sequence (recall that /, = /(m*) and 
that M* is the minimizer of / over all probality measures supported by S, see [H]). In 
addition, 

hm («„) = I{^lo) - CmaxMo(^) = , 

which combined with (|5.4ip and (|5.9p yields (|5.16p . 
Next it remains to prove that 

^n^CUrlj(Wn) Mo 

n — >-+oo 

strongly in {Cq'^{A)Y . Indeed, given an arbitrary function (p G Cq'^{A)^ we have 

J A J A J A 

-(^0,^) + — / (k„p-l)j«)-V^^ +0(1), 

J A 

where we have used (j5.13p and the fact that curlj'(w*) = in ^. Arguing as in the proof 
of Proposition 15. 21 we estimate 



1 



< c'-^\\j{uI)\\^^^a)^/eZM < = 0(1) . 



l^ii J A 

where the constant C only depends on ip, and the proof is complete. □ 
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